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Abstract 

Models for branched networks are often expressed as the minimization of an 
energy over vector measures concentrated on 1-dimensional rectifiable sets with 
a divergence constraint. We study a Modica-Mortola type approximation in¬ 
troduced by Edouard Oudet and Filippo Santambrogio, which is defined over 
vector measures. These energies induce some pseudo-distances between functions 
obtained through the minimization problem min{M“(M) : V • u = /'*" — /”}. We 
prove some uniform estimates on these pseudo-distances which allow us to establish 
a T-convergence result for these energies with a divergence constraint. 


1 Introduction 

Branched transportation is a classical problem in optimization: it is a variant of the 
Monge-Kantorovich optimal transportation theory in which the transport cost for a mass 
m per unit of length is not linear anymore but sub-additive. More precisely, the cost 
to transport a mass m on a length I is considered to be proportional to for some 
a €]0,1[. As a result, it is more efficient to transport two masses mi and m 2 together 
instead of transporting them separately. For this reason, an optimal pattern for this 
problem has a “graph structure” with branching points. Contrary to what happens in 
the Monge-Kantorovich model, in the setting of branched transportation, an optimal 
structure cannot be described only using a transport plan, giving the correspondence 
between origins and destinations, but we need a model which encodes all the trajectories 
of mass particles. 

Branched transportation theory is motivated by many structures that can be found 
in the nature: vessels, trees, river basins... Similarly, as a consequence of the economy 
of scale, large roads are proportionally cheaper than large ones and it follows that the 
road and train networks also present this structure. Surprisingly the theory has also had 
theoretical applications: recently, it has been used by F. Bethuel in jj] so as to study the 
density of smooth maps in Sobolev spaces between manifolds. 

Branched transportation theory was first introduced in the discrete framework by E. 
N. Gilbert in |14| as a generalization of the Steiner problem, fn this case an admissible 
structure is a weighted graph composed of oriented edges of length li on which some mass 
mi is flowing. The cost associated to it is then limf and it has to be minimized over 
all graphs which transport some given atomic measure to another one. More recently, the 
branched transportation problem was generalized to the continuous framework by Q. Xia 
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in |23| by means of a relaxation of the discrete energy (see also EH). Then, many other 
models and generalizations have been introduced (see m for a Lagrangian formulation, 
see also for different generalizations and regularity properties.). In this paper, 

we will concentrate on the model with a divergence constraint, due to Q. Xia. However, 
this is not restrictive since all these models have been proved to be equivalent (see jSj 
and [IS]). 

In this model, a transport path is represented as a vector measure u on some open 
set Q C such that V • tt = for two probability measures and . Then 

the energy of u is defined as M°^{u) = if n is a vector measure concentrated 

on a rectifiable 1-dimensional set M on which u has multiplicity 9 w.r.t. the Hausdorff 
measure (see |3| for more details). In this framework, u must be considered as the 
momentum (the mass 9 times the velocity) of a particle at some point. Then (V • u){x) 
represents the difference between incoming and outcoming mass at each point x. 

In this paper, we are interested in some approximation of branched transportation 
proposed by E. Oudet and F. Santambrogio few years ago in |T 8 | and which has interesting 
numerical applications. This model was inspired by the well known scalar phase transition 
model proposed by L. Modica and S. Mortola in |16| . Given u G i?^(H,M'^) for some 
bounded open subset H C E. Oudet and F. Santambrogio introduced the following 
energy: 

Mf(u)=e-T'i f \u\f^ + e^^ f \Vu\\ 

Jn Jn 

where /3 € (0,1) and 71,72 >0 are some exponents depending on a (see (|2.4p l. If u does 
not belong to the set the value of M“ is taken as -|-oo. 

We recall the heuristic which shows why M“ is an approximation of M“ (see my- 
assume that /i (resp. is a point source at Si (resp. S' 2 ) with mass m. Then, it 
is clear that the optimal path for M" between these two measures is the oriented edge 
S = {Si,S2) of length I with a mass m flowing on it. We would like to approximate 
this structure, seen as a vector measure u concentrated on S, by some vector fields v 
which are more or less optimal for M“. What we expect is that v looks like a convolution 
of u with a kernel p depending on e and m: v = u* pn for some R = i?(e, m), where 

Pr{x) = R-'^p{R-^x) (1.1) 

for some fixed smooth and compactly supported radial kernel p € (^^(M*^). Then the 
support of V is like a strip of width R around S so that |u| is of the order of m/R'^~^ and 
|Vu| is of the order of m/R^. This gives an estimate of M“(u) like 

A4“(u) ~ e-'^^R^-^{m/R'^-^)^l + e'^^R^-^{m/R'^)^L (1.2) 


With our choice for the exponents 71 , 72 and (3, the optimal choice for R is 


1—7 

R = , 


(1.3) 


where 


7 = 


72 

d + l' 


(1.4) 


2d-P{d-l) 

This finally leads to M“(u) ~ m" as expected. 

It was proved in |18| that, at least in two dimensions, the energy sequence (M")£>o 
T-converges to the branched transportation functional cqM® for some constant cq and 
for some suitable topology (see Theorem 12.11 page [5|). This result has been interest¬ 
ingly applied to produce a numerical method. However, rather than a T-convergence 
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result on M“ we would need to deal with the functionals M”, obtained by adding a 
divergence constraint: it should be shown that (u) := M^{u) + F-converges 

to cqM (u) := CoM^(u) + I'^.u=^+-^-1 where /g G is some suitable approximation 
of /i"*" — ii~ and Ia{u) is the indicator function in the sense of convex analysis that is 
0 whenever the condition is satisfied and +oo otherwise. Even if this property was not 
proved in |18| , the effectiveness of the numerical simulations made the authors think that 
it actually holds true. Note that an alternative using a penalization term was proposed 
in | 20 ) to overcome this difficulty. 

In section [2] we recall Xia’s formulation of branched transportation and its approx¬ 
imation M" introduced by E. Oudet and F. Santambrogio. The longest part of this 
paper, section [31 is devoted to a local estimate which gives a bound on the minimum 
value ■= min{M"(u) : V • u = /} depending on ||/||ii, ||/||l 2 and diam(n) 

(see Theorem 13.21 page (61). In section U we deduce a comparison between (i“ and the 
Wasserstein distance with an “error term” involving the norm of /'*' — f~ . As an ap¬ 
plication of this inequality, in the last section, we will prove the following T-convergence 
result which was lacking in |18| 

Theorem 1.1. Let {fe)e>o C L?‘{LL) be a sequence weakly converging to /i as measures 
when e ^ 0. Assume that the sequence {fe)e>0 satisfies 

[ fs{x)dx = 0 and e^ 2 ||j^|| 2 ^ —^ o_ 

Jn 

There exists a constant cq > 0 such that the functional sequence T-converges to 

cqM as e —>■ 0. Moreover cq is the minimum value for the minimizing problem dST]). 

This answers the Open question 1 in [201 [T 8 ] and validates their numerical method. 


2 Mathematical setting. 

The branched transportation energy In all what follows, we will use the model 
proposed by Q. Xia (see |23) and [24] i: 

Let d > 1 be an integer and be some open and bounded subset of M'^. Let us 
denote by Aidivi^) the set of finite vector measures on O such that their divergence is 
also a finite measure: 

■Mdivi^) ■= 1^^ measure on Ll valued in < -|-oo| , 

where := \u\{Q) |V • u|(Ll) with 

|n|(n) := sup tf ■ du : G C{Ll, M'^), ||V’||l°° < 1 

and, similarly, 

|V • u\{Ll) := sup • du : cp ^ C^(n,M), ||v3||l°o < 1 

In all what follows, V • u has to be thought in the weak sense, i.e. f (/^V - u = — f V(/J • dtt 
for all ip G C^(0). Since we do not ask (p to vanish at the boundary, V • u may contain 
possible parts on dLl which are equal to u • n when u is smooth, where n is the exter¬ 
nal unit normal vector to dTl. In other words, V • u is the weak divergence of u \ q , in 
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where 1 q is the classical indicator function of n, equal to 1 on and 0 elsewhere. 
From now on, the notation lx for the classical indicator function of a set X and Ix 
for the indicator function in the sense of convex analysis (equal to 1 inside and +oo 
outside) will be used. is endowed with the topology of weak convergence on 

u and on its divergence: i.e. Un ^ u if Un ^ u and X-Un X-u weakly as measures. 

Given 0 < a < 1, the energy of branched transportation can be represented as follows 
for measures u £ M-div 

/ N f if tt can be written as m = C/(M, 0 , 

M"(n) = ^ (2.1) 

1 +00 otherwise, 

where U{M,9,^) is the rectifiable vector measure u = 9^ ■ with density 9^ with 
respect to the Hausdorff measure on the rectifiable set M. The real multiplicity is 
a measurable function 9 : M ^ M'*' and the orientation ^ : M ^ S'^~^ C is such that 
^(x) is tangential to M for "H^-a.e. x E M. 

Given two probability measures /r'*' and on Q, the problem of branched trans¬ 
portation consists in minimizing M“ under the constraint X ■ u = /r'*' — /r“: 

inf{M"(u) : u £ Aidivi^) and V • u =. (2.2) 

Note that, if fj,^{dQ) = 0, the divergence constraint implies a Neumann condition on u: 
tt • n = 0 on 912. 


Functionals M“ For the minimum value in (|2.6I) to be finite whatever /i+ and 
in the set of probability measures, we will require a to be sufficiently close to 1. More 
precisely, we make the following assumption: 

l-i<a<l. (2.3) 

Q. Xia has shown in |23| that, under this assumption, there exists at least one vector 
measure u £ Aidivi^) such that M°‘{u) < -|-oo. 


We are interested in the following approximation of M“ which was introduced in m 
for all u £ Aidiv{^) and for all open subset w C 12, 





if tt E 
otherwise. 


where /3, 71 and 72 are three exponents depending on a and d through: 


2 — 2d + 2ad 
3 — d + a{d — 1) 


(2.4) 


and 

7 i = (d — 1)(1 — a) and 72 = 3 — d -|- a{d — 1). 

Note that inequality 1 — l/d<a<l implies that 0 < /? < 1. When a; = 12, we simply 
write 

M“(u,f2) =:M“(n). 
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(2.5) 


We point out the 2-dimensional case where M" rewrites as 

M"(ri) = f \u{x)\^ dx + f |Vtt(x)pdx, 

Ju Jo. 

where /3 = 

Given two densities /■*■, f~ G L^(0) := {/ G L?'{VL) : / > 0} such that f f~^ = J f~, 
we are interested in minimizing M^{u) under the constraint V • u = 

inf |M"(u) : u ^ H^{Q) and V-u = /^ —/“}. (2.6) 

The classical theory of calculus of variation shows that this infimum is actually a mini¬ 
mum. A natural question that arises is then to understand the limit behavior for mini- 
mizers of these problems when £ goes to 0. A classical tool to study this kind of problems 
is the theory of T-convergence which was introduced by De Giorgi in |12| . For the def¬ 
inition and main properties of F-convergence, we refer to m and [8]. In particular, if 
M“ F-converges to some energy functional Mq and if (ug) is a sequence of minimizers 
for M“ admitting a subsequence converging to u, then, tt is a minimizer for Mq. By 
construction of M“, we expect that, up to a subsequence, M" F-converges to cqM". In 
the two dimensional case, we have the following F-convergence theorem proved in |18) : 

Theorem 2.1. Assume that d = 2 and a G (1/2,1). Then, there exists a constant c > 0 
such that (M“)g>o T-converges to cM°^ in Mdiv{^) when £ goes to 0. 

Nevertheless, this does not imply the F-convergence of M^{u) -|- to 

M^(u) Indeed, the F-convergence is stable under the addition of con¬ 

tinuous functionals but not l.s.c. functionals. Gonsequently, we cannot deduce, from this 
theorem, the behavior of minimizers for (12.6p . For instance, it is not clear that there 
exists a recovery sequence (ue), i.e. Ug converges to u in Aidivi^) and M^{us) converges 
to M°‘{u) as e —)• 0, with prescribed divergence V • rtg = /■*■ — f~. To this aim, we require 
some estimates on these energies and this is the purpose of this paper. 

Distance of branched transportation We remind our hypothesis 1 — 1/d < a < 1. 
In [23], Q. Xia has remarked that, as in optimal transportation theory, M“ induces a 
distance d" on the space V{^}) of probability measures on fl: 

d"(^’'',/r“) = inf {M"(ri) : u £ Aidivi^) such that V • n =/r’*'—, 

for all G Thanks to our assumption a > 1 — 1/d, d" is finite for all 

G V{0.) and it induces a distance on the set P(fl) which metrizes the topology of 
weak convergence of measures. Actually, d“ has a stronger property which is a comparison 
with the Wasserstein distance: 

Proposition 2.2. Let //■*■ and be two probability measures on Q. We denote by Wp 
the Wasserstein distance associated to the cost {x,y) ^ |x — y\P for p > 1. Then, one 
has 

Vki/„(/i+,/r-) < d“(;U+,/r-) < C 

for a constant C > 0 only depending on d, a and the diameter of LI. 

We refer to m for a proof of this property (see also |3|, and ^ for an alternative proof) 
and [22], [21] for the definition and main properties of the Wasserstein distance. In the 
same way, we define d“ as follows: 

d“(/+,/“) = inf {m“(u) : u G dd^IE^'^) such that V-u = /+-/'}, (2.7) 
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where f^,f G -^+(^) satisfy /q ■ Although d°‘ is a distance, it is not 

the case for d" which does not satisfy the triangular inequality. Actually, because of the 
second term involving |Vttp, M" is not subadditive. However, for uj,... ,Un in A4div(^), 
the inequality |Vtti + • • • + < n{|Vuip + • • • + |VunP} implies 

( n \ n 

) < n^M"(ui). 

i=l J i=l 

In particular, d" is a pseudo-distance in the sense that the three properties in the following 
proposition are satisfied: 

Proposition 2.3. Let /+, f~ and fi,..., /„ be densities, i.e. non negative 
funetions whose integral is equal to 1. Then one has 

1- = 0 implies /+ = 

dnf+j-) = dnf-,n, 

3 . d“(/o, fn) < n [d“(/o, /l) + d“(/i, /2) + • • • + dfifn-lJn)] • 

3 Local estimate 

We remind our assumption (12.31) which insures that d°‘{^i^, fi~) is always hnite. Our goal 
is to prove that d" enjoys a property similar to the following one. 

Proposition 3.1. Let Qq = (0, C 6 e a cube of side length L > 0. There exists 
some constant C > 0 only depending on d and a such that for all non negative Borel 
finite measure p, of total mass d > 0, 

d^{p,0ho) < C rL, 

where do is the Dirac measure at the point cq^, the center ofQo. 

Since d“ (/■*■, /“) is only dehned on functions /^, to do so, we first have to replace 
96q by some kernel which concentrates at the origin when e goes to 0. Let p G 
be a radial non negative function such that p = 1, where B C is the unit ball 
centered at the origin, and define pg^s '■= PR in m, where 

R=: Rg^^ = e'^e^. 

Here, we recall that R and 7 = were introduced in dn. Let Q be a cube in 
centered at some point cq G and / G L^(Q) be a density such that Jg / =: 6q. 
Then, we will denote by pq the kernel Opg^e refocused at cq with a small abuse of notation 
(indeed, pq also depends on /): 

Pq{x) = OqpoQ^six - cq). 

The main result of this section is the following theorem 

Theorem 3.2 (Local estimate). Let us set Qo = (0,L)'^ for some L > 0. There exists 
C > 0 only depending on a, p and d such that for all f G L‘^{Qq) with f = 6, we 
have 


6 


If supp Pqq C Qo then, there exists u E Hq{Qq) such that ■ u = f — pq^ and 
dfif,pQ,)<M^iu)<C{e^L + e'^^\\f\\l,} and \\uhi<CLe. 


• Otherwise, there exists u € Hq{Qq) such that 

d<f{f,pQ,)<M^{u)<Cs'^^f\\l, and ||nbi<CL0, 
where Qo = 2supp/9Qo := B{cQ^,2Re^s). 

Remark 3.3. The Dirichlet term, ||/||^ 2 ) in the estimates above is easily understand¬ 
able. Indeed, if e is very large so that one can get rid of the first term in the energy M“, 
then, one can use a classical Dirichlet type estimate, that is Theorem 13.41 below. On the 
contrary, for e very small, the T-limit result on M“ tells us that these energies are close 
to M" so that it is natural to hope a similar estimate as the one for M": that is to say 
an estimate from above by (see El)- 

The main difficulty to prove Theorem 13.21 is the non subadditivity of the pseudo¬ 
distances A“. Indeed, our proof is based on a dyadic construction used by Q. Xia in 
|23| to prove Proposition 13.11 (see also m)- This gives a singular vector measure u which 
is concentrated on a graph. Since M" contains a term involving the norm of Vtt, 
we have to regularize u by taking a convolution with the kernel pe^e on each branch of 
the graph {9 being the mass traveling on it). Unfortunately in this way, two different 
branches are no longer disjoints. 

It is useful to see that we have a first candidate for the minimization problem (12.71) . 
This candidate is of the form u = V()), where (f is the solution of the Dirichlet problem 

lAi, = /+ - /- in Q. 

I 0 = 0 on dQ. 


Then n = V0 satisfies V-u = in Q and u{x) E Mn a.e. on dQ where n stands for 

the external unit normal vector to dQ. Alternatively, one could consider Neumann ho¬ 
mogeneous boundary conditions for cf rather than Dirichlet boundary conditions. Then, 
one would obtain u{x) • n = 0 a.e. on dQ. Theorem 13.41 below gives a better result in the 
sense that the candidate u vanishes at the boundary: 

Theorem 3.4. Let Qo = (0, L)'^ he a cube of side length L > 0. There exists C > 0 
only depending on d such that for all f E Lq{Qo), there exists u E 77q(Qo)1^^) solving 
V ■u = f and satisfying ||n||ii(QQ) < CL ||/||l 1 (Qo) together with 

I |Vn|M <C||/|U2(Q„), 

Qo / 

where LKQo) = ^f £ L'^iQo) : /(x) dx = o|. 

For a proof of this theorem, see, for instance, Theorem 2 in [7]: the only difference with 
Theorem 13.41 is that we add the estimate ||ti||Li(Qg) < CL ||/||l 1 (Qq) which can be easily 
obtained following the proof of J. Bourgain and H. Brezis. The corresponding property 
formulated on a Lipschitz bounded connected domain D is also true (see Theorem 2’ in 
m) except that the constant C could depend on D in this case. 

Of course, this candidate is usually not optimal for (12.71) and this does not allow for 
a good estimate because of the first term in the definition of M". For this reason, we 
have to use the dyadic construction of Q. Xia up to a certain level (“diffusion level”) from 
which we simply use Theorem 13.41 
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3.1 Dyadic decomposition of Qq and “diffusion level” associated to f 

Let us call “dyadic descent” of Qo = (0,L)'^ the set Q = Uj>o where Qj is the 
“dyadic generation”: 

Qj = {{xi,... ,Xd) + 2~^Qo : Xi & {k2~^ L : 0 < k < 2^ — 1} for i = 

Note that Card(Qj) = 2^^. For each Q ^ Q, let us define 

• 'D{Q): the descent of Q, the family of all dyadic cubes contained in Q. 

• A{Q): the ancestry of Q, the family of all dyadic cubes containing Q. 

• C{Q): the family of children of Q composed of the 2'^ biggest dyadic cubes strictly 
included in Q. 

• F{Q): the father of Q, the smallest dyadic cube strictly containing Q. 

We now remind the dyadic construction described in |23) which irrigates / from a point 
source. We first introduce some notations: fix a function / € L^(Qo) with integral 6 
and let <5 E Q be a dyadic cube centered at cq E Then we introduce Oq the mass 
associated to the cube Q as 



If 9q a 0, we also define the kernel associated to Q through 

Vq{x) = Pr{x), (3.2) 

where is defined in (II.ip for 

R = RQ-.= e-’sf' , 7 =^ . 

Here 7 was dehned in Define also the weighted recentered kernel by 

Pq{x) = 9qPq{x - cq) (3.3) 

if 7 ^ 0 and Pq{x) = 0 otherwise. Lastly, we introduce the point source associated to 
the cube Q as 

Sq := Oq X Dirac measure at point cq. 

We are now able to construct a vector measure X such that M“(X) < + 00 . First define 
the measures Xq as below: 

^Q= ^ ^Q'^Q''^\cq,cq,P (3-4) 

<3'eC(Q) 

C-Qf — Cq 

where nQ/ = - 7 -. Then, we have 

\\cq'-cq\\ 

^ ' ^Q= Sq' - Sq 

Q'eCiQ) 

and the energy estimate 

M"(Xq) < 2'^-20“diam(Q), 

where diam((5) stands for the diameter of Q. Finally, the measure X = YIqgq^Q solves 
X ■ X = f — Sqq and satisfies 

M“(X) < C'0“diam(Qo)- 


Indeed, it is enough to apply the following lemma with A = a: 




Lemma 3.5. Let Q G Q and A e]1 — 1/d, 1]. There exists a constant C = C{\,d) such 
that 

0Q/diam((50 < diam((5). 

Q>&V{Q) 

Proof. Let jo > 0 be such that Q G The definition of P(Q), the Jensen inequality 
and the fact that d — 1 — Ad < 0 give 

9^, diam(Q') = Y diam(Q) Y ^Q' 

Q'eViQ) j>o Q'eViQ)nQj^+j 

< diam((5) E 2-j2^d 2-Ti Y ^Q' 

i>o \ Q'6i’(Q)nQjo+j 

< 9 ^ diam(g)y^2^('='-i-^'^) 

j>0 

< Cdg diam((5). □ 

Now, the idea is to replace each term in (13.4p by its convolution with the kernel pg'- 
Unfortunately, this will make appear extra divergence terms around each node. We have 
to modify X so as to make this extra divergence vanish using, for instance, Theorem 13.41 
Furthermore, we cannot follow the construction for all generations j > 1, otherwise the 
“enlarged edges” (convolution of a dyadic edge and the kernel pe.e) niay overlap. This is 
the reason why we introduce the following definition: 

Definition 3.6 (“Diffusion level”). For a cube Qq and / G L^(Qo) we define the set 
VlQo, f) or 'D(f) C Q as the maximal element for the inclusion in the set 

A = {D C Q : \/Q G D, A{Q) UC{F{Q)) C D and supppg C Q} . 

If A = 0, that is supp pg^ ^ Qq, we take the convention iD(/) = 0. For all x G Qq, define 
also the “generation index” of x associated to / as 

j(/,x)=max{j : 3Q G P(/) n Qj, x G Q} G N U {±oo}, 

where the convention max(0) = —oo has been used. 

In this way, each cube in T’(/) contains the support of its associated kernel. Moreover, 
if Q is an element of T>(/), then all its ancestry and its brothers (i.e. elements of the 
set C{F[Q))) are elements of T>(/). T’(/) can be constructed by induction as follows: 
first take j = 0 and 'D(/) = 0. If supp pg^ C Qq then add Qq to the set T’(/) and j is 
replaced by j + 1. For all cubes Q in A n Qj-i'- if all cubes Q' G C{Q) C Qj are such 
that their associated kernels are supported on Q' then T’(/) is replaced by 'D{f) DC{Q). 
If T’(/) has been changed at this stage j is replaced by j + 1 and the preceding step is 
reiterated. This process is repeated for j > 1 until it fails. 

Let T>min{f) be the set of all cubes in T’(/) which are minimal for the inclusion. If 
Fminif) / 0, we also define 

D{f)= U Q ■ 

Note that this is actually a disjoint union: two distinct cubes in 'Dmin{f) are disjoint. 
Indeed, for Q, Q' G 'Dminif) C Q, either QdQ' = ^ or Q and Q' are comparable: Q C Q' 
or Q' C Q. In the last case, since Q and Q' are minimal, we deduce that Q = Q'. 
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Moreover, it is not difficult to see that, if Pmm(/) 7^ 0) then D{f) = {x £ Qo : 
j{f,x) is finite} and also that /(x) = 0 whenever j{f,x) = +oo, where / is the precise 
representative of / (i.e. the limit of the mean values of / on small cubes). Indeed, assume 
that Q G T^if) is a cube of side length Lq. Then, by definition, supppg C Q and for 
some constant C depending on p and for v = one has < CLq and so 

Since 1/u — d = jg positive, we deduce that Lq cannot be arbitrarily small 

if there exists x £ Q such that f{x) > 0. Moreover, if f{x) > r] a.e. for some p > 0, then 
there exists some constant Crj > 0 depending on p, e, d and a such that 

VQ G P(/), Lq > Cr,. (3.5) 

In particular, one can deduce that T>min{f) = 0 if and only if T>{f) = 0 or f{x) = 0 
a.e. Indeed, if T’(/) = 0, then it is clear that Prnm(/) = 0- Conversely, assume that 
L’minif) / 0 (he. Qo £ T>(/)) and that there exists x £ Qq such that f{x) > 0. Since 
UgeCl/) negligible for the Lebesgue measure, one can assume that x £ UQeD(/) Q- 
Then 0 < j{f,x) < +oo and so there exists a minimal cube Q £ T>{f) containing x. 
Then Q £ 'Dminif)- Indeed, if Q' £ T>{f) and Q' ^ Q, then A{Q) C T’(/) and there 
exists Q" £ A{Q) such that Q" C Q and x £ Q" which is a contradiction. 

We are now able to define two approximations of / which are useful for our problem. 
The first is a dyadic approximation of / by an atomic measure, 

( Y. if^min(/)/0, 

Ae/ = < Q&Vminif) 

I Sqo otherwise, 

where we recall the definition of Sq := 9q5cq- We also define an approximation in 

hHQo), 

I Y '^min(f) A 0, 

{ PQo otherwise, 

where pQ is defined in (13.311 . The following result shows in which sense X^f is an approx¬ 
imation of / and justifies the term “diffusion level”. Indeed, this proposition indicates 
that we get a good estimate by using a local diffusion from X^f to /, i.e. minimizing 
Jg |Vnp over the constraint V • n = X^f — f for all Q £ 'Dminif) (see Theorem EaD. 

Proposition 3.7. There exists a constant C > 0 depending on d and p such that for all 

f G LliQo), 

d“(A,/,/) + d“(A,/,/) <Ce^^f\\YQ^y 

More precisely, if supp pQ^ C Qo, there exists u £ Hq{Qo) such that V ■ u = f — X^f as 
well as 

Mf‘{u) < Ce'^^\\f\\l 2 and ||n||ii < C diam(Qo)||/||Li- 

If supp pQg ^ Qo the same estimates hold but the condition u £ Hq{Qo) has to be replaced 
by u £ Hq{Qo), where Qo is a cube containing Qo and supp/og^. 

Proof. First assume that supp/Og^ C Qo i.e. Qo £ 'L>{f). If Vminif) = 0) then f{x) = 0 
a.e. and the proposition is trivial. Hence, one can assume that 'Dminif) 7^ 0- Then / 
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is supported on D{f) and 'Dminif) =• {Qi}iel is a finite or countable partition of D(f). 
Denote for simplicity Di := diam((5j), fi ■= (restriction of / to Qi), Oi := 6Q^ and 
Pi ■= PQi= PRi for ^ where 

R, := Rq^ = 

Since Qi is minimal in 'D(/), we deduce that, for some constants C, C > 0, 

C'Ri <Di< CRi. (3.6) 

Indeed, the first inequality follows from the fact that suppp* C Qi and diam(supppj) = 
cRi for some constant c depending on p. For the second inequality observe that, since Qi 
is minimal, there exists Q G C(Qi) such that supp pq^Q and hence Rq > d diam((2) = 
dl2Di for some constant c' > 0 depending on p. Since 9q < 9q^ = 9i, one has Rq < Ri 
and the second inequality follows. 

Now, Theorem 13.41 allows us to find Ui G HQ{Qi) such that V ■ Ui = pi, < 

C\\gi\\L^{Qi) ^rid \\ui\\Li(Qi) < C DiWgiW^^Q.p where pi := fi - pi. Since m vanishes 
at dQi, one can extend Ui by 0 out of Qi to get a function in Ff^(M‘^): for the sake 
of simplicity, this function is still denoted by Uj. Consequently, u = 'Yhi'^i belongs to 
Hq{Qq) and S/ ■ u = f — Xef ■ It remains to estimate M^{u) and ||rr||_Li(( 3 Q). First of all, 

||w||i,i(Qo) < Wuih^iQi) < C diam(go) Y WdiW^m 

i i 

and the inequality ||5i||Li(Qi) ^ leads to ||rt||x,i < 2C diam(Qo)||/||Li ^ required. 

Let us compute the L^-norm of pp 

llftllb(Q.) = Oi ll/’R,llb(Q.) = oIrPMIhq.) = (OOlRp- 

By a Cauchy-Schwarz inequality, 

/.)' < mmhm = ( 3 . 7 ) 


which, together with (13.6p . gives 

Mhm < cRfm^Q^^R-^ = c ml^^Q^y 

Since ||ui||j^i(Q.) < CH/i - we get ||ni||j^i(Q.) < C'||/i||i 2 (Q,). Now, because the 

energy is local and since each Ui is supported on one has 


n ^ / /* r 

A4“(n) = M“(u,) = / In/ + / | Vu,| 

i=i i=i V 

By construction of Ui, one has 

[ IVUjp < ll'Will^qQ.) < ^llfl'*lli2(Qi) — (ll7'*lli2((3i) + ll/*lli2(Qi)) — 

j Qi 

It remains to estimate the first term. First of all, we use the Holder and 
inequalities as follows: 

[ !«/ < \Qp-oi'^ 

Qi 


/3/2 


<Di 


d-d/3/2 


Di 


iVud 


/3/2 


< DWi\ 


2 

D{Qi) ■ 
Poincare 


P 

D{QiV 
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where v = /3 + d — In view of (13.7p and (13.6p , we have 


1 —y 
) d-i 


1-7 


A < CR, = Ce^dt^ < WMlHQ.))^ 


and, introducing 5 := 1 — 




Finally, since -71 + ^ = 72 and /3 + = 2, we get 


v(^ —r) 


--71 




[ \Ui\^ 

Qi 

The proof of the second inequality is quite similar but easier: 

n n 

d“(A./, /) < E /*) < OfR 


i=l 


i=l 


/dlL2(Q.). 


(3.8) 


Once again, applying (13.7D and then (|3.8I1 . we get 

d“(AA,/)<CA^||/||i.. 

1-7 

In the case where supppQp ^ Qo; he. Rq^ := > CL [L being the side length 

of Qo and C a constant depending on p), the proof is the same. Indeed, we just apply 
Theorem 13.41 to g = f — pg^ and the same computations as above lead to the same 
result. □ 


3.2 Proof of Theorem 13.21 

Let Qo = (0,T)'^, L > 0 and / G L‘^{Qq) with f = 0. In the case where supppQg ^ 
Qo, Theorem 13.21 is a particular case of Proposition 13.71 Consequently, one can assume 
that supppQg C Qo i.e. Qo G Rif)- In the case where V[f) = {Qo}) one has Ag/ = 
and Theorem 13.21 is a consequence of Proposition 13.71 as well. For this reason, one can 
assume that C{Qo) C R{f)- Moreover, up to replacing / by f + p for some small constant 
p > 0 and passing to the limit when p —)• 0, one can assume that 'P(/) is finite. Indeed, 
in view of (13.511 . 'D{f + p) is hnite since for all Q G 'D{f + p), diam(Q) > > 0. 

Our aim is to prove that there exists C > 0 only depending on a, d and p such that 

The idea of the proof is to approximate the vector held X = '^ Xq of the previous 
section (see (13.40 by a vector held in using the kernel p. In this part, we will use the 
notations of the previous section: in particular, the dehnition of T>{f) in Dehnition 13.61 
the measures Xq in (13.4p and X = Ylgevif) ^Q- 

Since C(Qo) C Rif), we can construct the regularized vector held Y by the formula 

E 2e. 

Q&nf) 

Q^Qo 

where, for all Q G R{f) such that Q 7^ Qo (see Figure [1]), 

Zq := 0Q ng Pq * (3-9) 
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nq being the normalized vector nq 


\\cl-l7Z\\ denned in O. 



Figure 1: A square Q and its 4 dyadic children Qi with the associated vector field Zq 

By definition of the kernel 'pq, one has 

M^{Zq) < Ce^ diam(g). (3.10) 

This a consequence of the choice of Rq as a minimizer in (ll.2p . Indeed, for the sake of 
simplicity, let us assume that supp p is the unit ball centered at the origin. Then Zq is 
concentrated on a strip of width Rq around the segment S = [cF{q)jCq], i.e. 

supp^Q C S := {x G : dist(x,S') < Rq} (3-11) 

and Zq satisfies the two estimates 

\\Zq\\L^ < CeqRl^-'^ and W^ZqU^o < COqRQ^ (3.12) 

Then, the same computations as in (11.211 and the fact that Rq < diam(Q) give (13.1011 . 
Let us estimate the norm of Y which has to be controlled by 9 as stated in Theorem 

[T2) 

E l|2«lltHS)SE E «0t2->=L9. 

i>i Q&v{f)nQj j>i QeVif)nQj 

Note that 

V-Y = pq,-h- XJ, 

where h stands for the extra divergence, h can be written as 










where Pq'^q represents the kernel pQ/ translated at cq, center of Q, and, for the sake of 
simplicity, the set of all cubes Q such that C{Q) C 'D{f) has been denoted by 'Dfr,- 

Vfr. := {Q G V{f) : C{Q) C P(/)}. 

Since V • y = pq^ — f + {f — Xef) — h ^ pq^ — /, we have to slightly modify the vector 
field Y. This will be done replacing Y by 

y = y + Fi + Fa, 

where Vi, V 2 G are constructed so that V ■ Vi = h and V • V2 = Ag/ — /. 

The construction of Vi and the estimate of M“(Vi), ||Vi||2,i will be the object of the first 
step. In the second step we prove that M“(y) < CO^L. Then, Proposition 13.71 allows 
us to construct Va G such that V • Va = A^/ — / with an estimate on M"(Va) and 

r2||Ll. 

First step: Correction at the nodes, constrnction of Vi. For all Q G 'Djr., let 
Bq be the support of pq. Since suppp has been supposed to be the unit ball centered 
at the origin and pq{x) = 6qpjig{x — cq), we have Bq = B{cq,Rq) C Q. Let us dehne 
the extra divergence corresponding to this node, 

hq= pq- ^ Pq',q ■ 

Q'eC(Q) 

For each Q G P/,.., thanks to Theorem 13.41 we can hnd Vq G H^^Bq) such that V • Vq = 
hq and HVQllHbBQ) ^ C \\^q\\l^{Bq)- this case, because hq is radial up to a 

translation, we essentialy use the proposition in dimension 1 which is quite easy and 
gives better estimates. Let us give more details on this point: 

Lemma 3.8. Let d> 1 and B = B{0,R) C be a ball centered at the origin. There 
exists a constant C > 0 only depending on d such that the following holds: 

Let F G L°°(B) be a radial function: i.e. for a.e. x G B, F(x) = /(|x|) for some f G 
L°°(0,R). Assume that f^F = 0. Then, there exists a radial function V G 
such that V • y = F and 


I|vi^IIl°°(qo) - ^ II-^IIl°°(qo) • 

Proof. First of all, by a scaling argument, one can assume that R = 1. The vector held 
y : B dehned by y(x) = v{\x\)x for some Lipschitz function v : M'*' —>• M satishes 

V • y(x) = r^~'^[r'^v{r)]' 

in the distributional sense. Thus, if v is chosen as 

v{r)=r-‘^ r f{s)s‘^-^ds, 

Jo 

then y solves the following problem: 

f V • y(x) = F{x) on B, 

^ y(x) = 0 on dB. 

Moreover, for a.e. x G B, we have Vy(x) = v'{\x\)x (8) ify + u(|x|)Id, where Id is the 
matrix identity. In particular, we get ||Vy|| LjOO < C(||ry(r)|| ^00 + The second 

term in the RHS on the preceding equation is estimated by ||u||2,oo < = 

ll/lli^oo. For the hrst term, one has vfr) = —ds + r“^/(r) and so 
\\rv'{r)\\L^ < C||/||loc. Thus, ||Vy||Loc < C\\F\\l^. □ 
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Applying Lemma 1^81 to F = hg and R = Rq gives Vq G Hq{Bq) such that V • Vg = 
hq and 

< COqR-q^ , WVVqWl^Bq) < \Bq\\\VVq\\l^^Bq^ < COq. (3.13) 

Moreover, since Vg is supported on Bq = B{cq,Rq), we deduce that ||Vg||/,oo( 5 ^) < 
^QllVFgll ^ C'0gi?Q SO that Vg satisfies the same estimate as (|3.12p . In par¬ 
ticular, we get M“(Vg) < diam((5). Now define 

^1= E 

QeVfr. 

Since ||Vg||Li(Sg) < CRq\\Wq\\Bi(BQ) < c diam{Q)eq, Lemma [33] implies 

||lA||Li(go) < C'diam(Qo)6'<3o ^ C'^-^II/IUlQo) 

as required. Then, using the definition of M“ in (12.4p and the subadditivity of x ^ |x|^, 
one gets 

Mf(Vi)<e-^i ^ f \Vq\f^ + 2e^^ f ^ [V^g. iV^g], (3.14) 

Q&Vfr. ■’ ■’ Q,Q'&Vf^, : Q'cQ 

where A : B stands for the euclidian product of two matrices A = (Ajj)i<j B = 

{Bij ) i<i,j<d of size dx d: A : B := AijBij. For the estimate of |VI4ip, we have used 

the identity |VI/ip = VIA : VIA = 'Eq,q'&Vfath VVg : VVg/. Since Vg is supported on 
Q, VVg : VVg' vanishes except when Q n Q' / 0, i.e. Q C Q' or Q' C Q, thus justifying 
the factor 2 and the inclusion Q' G Q in (|3.14D . 

We need to estimate the two terms in (|3.14p . Since M^{Vq) < COq diam((5), thanks 
to Lemma 13.51 this term is less or equal than C6°‘L as required. Using the inequality 
WfgWF < ll/ll one can estimate the second term of (13.1411 by 

2e^^ ^ l|VUg|A^(BQ)||VUB^,|Ai(B^,). 

Q,Q'&Vfr. : Q'cQ 

Note that it would be more natural to use a Cauchy-Schwarz inequality {Lp‘-L?‘) at this 
step but, using it, we were not able to deduce the estimate by 6°'L. Once again, since 
Rq' < diam((5')) have 

llVUg'IAqs^,) < COq, < diwn{Q')R-}eq, = Cdiam(Q')e"^0Q7^. (3.15) 

Since 1 — ^ < 1 — ^5^ < 1, Lemma [33] gives 

E llV^Q'llLABg') < C£-^diam(g)0Q"^. 

Q'CVfr. ■■ Q'cQ 

Now, elementary computations on exponents a, 70 , 7 and Lemma [33] give successively 
7. = (d+l)7,o = 2-(d+l)^and 

E diam(Q)0gi2Q'^e-T'6(Q"^ = E diam(Q)6'§ < C0“L. 

Q^T)fr, Q^T^fr. 

Finally, we have obtained the desired inequality: M“(Vi) < C O^L. 
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Second step: estimate of the energy of Y on the node set. In order to get 
estimates on Y, it is convenient to divide Qq into 2 domains: the node set N and its 
complementary where 

N := IJ B{cq,cRq) 

Q&nf) 

and c > 0 is a constant which will be chosen later. By analogy with Vi, one can write 
Y|jv as a sum of vector fields Yq, where 

— i ^s(cQ,ci?Q) ~ if Q € '^fr. (see ([32])), 

1 '^b{cq,cRq)Zq otherwise. 

Now, from (I3.12p . we deduce the estimates (I3.13P satisfied by Vq are also true for Yq 
and consequently, we obtain M^(Y,N) < C 6°‘L as well (see (|2.4D for the definition of 
M^{Y,N)). 

Third step: estimate of the energy of Y out of the node set. Reminding that 

r= E ' 

Q^Qo 

considering that M“ is not subadditive (due to the term |Vyp), the first thing to do is 
to understand to which extent the supports of Zq can intersect. To this aim, let us note 
that if the constant c > 0 in (|3.2p is chosen equal to \/d or more, due to (|3.1ip . then each 
Zq restricted to iV'^ is supported on Q (see figured]): supp Zq n C Q. In particular, 
this implies that 

supp Zq n supp Zqi n N'^ 7 ^ 0 Q CiQ' ^ Q C Q' or Q' c Q. 

For this reason, M^{Y,N^) can be estimated exactly in the same way as we did for the 
estimate of M"(I/i) in (13.141) . Moreover, the Young inequality, ||/*/i||Li < ||/||£,i|/x|(M'^), 
valid for all / E L^(M'^), ^ E A4(E'^), and the definition of Zq in (13.9p . easily give 

llV^Q'lliRQ') < CeQ>R~} diam.{Q'). 

Since this estimate (which is the same as (|3.151) 1 and (I3.10p are the only ones we have 
used in the first step for the estimate of M“(Yi), we get M^{Y,N^) < C6°‘L as well. 

End of the proof of Theorem 13.21 Finally, the vector field V = Y + Vi + V 2 , where 
V 2 is given by Proposition 13.71 satisfies V • P = pQ^ — /, 

Mf(P) < 3{M“(y) + M“(Pi) + M“(P2)} < C{e<^L + e^^\\f\\l,} 

and 

IIpIIli < ||y||Li + llVhlUi + II^2 ||li < CLWfUi. 

4 Estimate between and the Wasserstein distance 

Our aim is to prove an estimate on the pseudo-distances similar to Proposition 12.21 
Because of the Dirichlet term in the definition of M", d" cannot be estimated only by 
the Wasserstein distance Wi but one has to add a term involving ||/+ — /~||i 2 . Using 
Theorem 13.21 we are going to prove the following theorem: 
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Theorem 4.1. Let Q = (0, L)'^ be a a cube of side length L > 0 in and e E (0,1). 
There exists C > 0 only depending on a, d and L such that for all f~^, f~ E L^(Q) with 
Jg /■*“ = Jq f~ = 1, there exists u E compactly supported on the set := {x E 

: dist(x, Q) < Ce'^} satisfying V ■ u = f := /'*' — f~ as well as 


dfif+,f-)<M^{u)<CH{wl-^^^-'^\f+,f-) + e'^^\\f\\l,) and ||n||^i < C, (4.1) 

where H : M'*' — > M"*" is the scalar function defined by H{x) = x + x^ for some A E (0,1) 
depending on a, and Wi stands for the Wasserstein distance associated to the Monge cost 
{x,y) |x - y\. 

Remark 4.2. One can replace the condition J /^ = 1 by f f^ = 9 > 0. Then, the 
constant C will also depend on 0: C = C{9,a,d,L). However, we can easily check that 
C is locally bounded with respect to 6, i.e. it is uniform for bounded values of 6. 

Remark 4.3. It is tempting to think that estimate dn]) also holds when H (x) = x which 
would be the natural choice. Indeed, if e is taken very small, since M“ F-converge to 
M" and because of Proposition 12.21 one can expect that df{f^,f~) ~ d'^{f~^, f~) < 
CWi{f~^,f~fi~^^^~°‘\ On the contrary, when e is very large, because of Theorem 13.4[ 
one can expect that d"(/+,/' -) c. e'T^fWl,. However, for technical reasons, due to the 
lack of subadditivity of the second term (Dirichlet energy) in the definition of M“, we 
were not able to reach the case H{x) = x. 

Proof. Our method to prove this proposition is an adaptation of that of J.-M. Morel and 
F. Santambrogio in m (see also Proposition 6.16. page 64 in [3]). 

Flp to replacing (/^,/~) by (/+ — /^ A f~,f~ — A /“), one can assume that 
/■*■ A /“ = 0, where for all x E Q, (/“ A f~^){x) = mf{f~{x), f~^{x)). Indeed, it is 
sufficient to note that, if are two measures with the same mass and z/ is a positive 
measure on Q then we have Wi{p~^ + u, p~ + v) = Wi{p^, p~). 

For the sake of simplicity, in all the proof, C > 0 will denote some constant only 
depending on a, d and L and big enough so that all the inequalities below are satisfied. 

Let /^, f~ E L‘^{Q) be two densities on the cube Q = (0, L)'^ such that Jq = 1. 
Chose an optimal transport plan H between /■*■ and f~ for the Monge-Kantorovich 
problem associated to the cost c(x, y) = \x — y\. Hence H satisfies the constraint P^H = 
f^{x) dx where P"*" (resp. P~) is the projection on the first variable x (resp. the second 
variable y) and dx is the Lebesgue measure. Moreover we have 


/ 

JQ 


\x-y\dU{x,y) = Wfif+,r)=:W. 


(4.2) 


So as to use the local estimate of the previous part, let us classify the set of ordered pairs 
{x,y) with respect to the distance \x — y\. More precisely, for j > 0, set 

Xj = {{x,y) £ : dj < \x - y\ < dj+i}, 

where dj = (2-^ — 1) re and w E (0,1) will 
is empty if dj > diam((5), i.e. j > J := 

restrict to integers j < J < (7(1 + | Intel): 

(14.2p immediately gives the estimate 

'^djej<W, where ej=U{Xj). (4.3) 

j 


be chosen later. In particular, do = 0 and Xj 
1112 ^ diam(Q) reason, one can 

we will assume that dj < diam((5). Moreover, 
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Next, for each integer j G consider a uniform partition of Q into cubes Qjk, 

k = 1,..., Kj, with side length dj+i- It is easy to estimate Kj by 

K, < CdJ^,. (4.4) 


For j > 0, set 


n, = n|x, ; Oj = n{x,y, ff = p±n, and /, = /+ - /r, 

Clearly, one has 

n = J^n, and /* = EC- 

j j 

In the same way, for j > 0 and 1 < k < Kj, set 

Iljfc = II|Xj n (QjfexQ) ’ ~ IIjfc(Q ) arid fjf. = Iljfc 

so that 

% = E "it; "j = E Of C = E ff- 

k k k 

lijk represents the part of the transport plan 11 corresponding to points in Qjk which are 
sent at a distance comparable to dj+i- In particular, is supported on Qjk and fjj. is 
supported on the cube Qjk with the same center but twice the side length of Qjk- As we 
did in (|3.3I1 . let us define pjk the kernel associated to Qjk by 

. , f p{Rjk{x - Cjk)) if 9jk / 0, 

= I 0 otherwise, 

1- 1 

where p G (M*^,M'*'), Rjk = P^9jk^ and Cjk is the center of Qjk- For the sake of 
simplicity, let us assume that supp p is the unit ball centered at the origin. Let Bjk := 
B{cjk,rjk) be the smallest ball containing Qjk and supppjfc = B{cjk, Rjk)- be. rjk = 
max{i?jfc, diam((5jfc)}. Thanks to Theorem 13.2[ it is possible to find a vector field Ujk G 
H^{Bjk) satisfying V • Ujk = fjk ■= fjk - fjk > \Wjk\\Lf(Bjk) ^ and 

M,k ■-= M^iujk) < C{ey,dj+, + e'ry\fjk\\l,^s^J. (4.5) 

Moreover, if Rjk > dj+i/2, the first term in the right-hand side of (14.5p can be omitted 
since one has 

9%dj+^<Ce^y\f,k\\h- (4.6) 

Indeed, in this case, writing 9 := 9jk and R := Rjk, one has 9‘^dj^i < 29^R and, using 
2 — a = we get 9°‘R = [9°‘~‘^R^~^'^\[9‘^R~‘^] = e'f^R~‘^9‘^. Then, (|4.6I) follows 

from the fact that, by the Cauchy-Schwarz inequality, we have 

j^-dg2 < f ^ 

k k 

Now, let us define the vector field u = '^^Ujk, which satisfies 

j,k 

y .u = Y^v ■ujk = Y,fjk = f ■= - r- 

j,k j,k 
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First note that 


u\\lI{Q) ^ — 2C'- 


In order estimate the energy of u, a similar development of (13.1411 and 

the Canchy-Schwarz ineqnality give 

J /Kj 

i=i \k=i 


E ^/WkVWi 

k {k,l)&Ij 


(4.7) 


where Ij stands for the set of pairs {k, 1) satisfying k ^ I, 6jk > Oji and Bjk H Bji ^ 0. 
We have to estimate the two terms in the right-hand side of (|4.7F 


Estimate of the first term in (14.7p We recall that Mjk < 6‘^f^djj^i + 

For the second term, note that 

(4-8) 

Indeed, since A f~ = 0, for all j, k, one has A fjj, = 0 as well. In particular, 

, ll/lli 2 (Q) = II/’'"IIl 2 (q) + 11/ lli 2 (Q) and ([4^ 
follows from the super-additivity of the power function x —>■ \x\^ for p > 1 : \x + y^’ > 
|x|^ + \y\^ for X, y G M whenever xy > 0. 

For the first term, applying successively the Jensen inequality with power a. G (0,1), 
the Holder inequality, (14. 3 p and the fact that Kjdj+i = Cdj^f (see (I4.4p h one gets 

< J2dj+iKj[0j/Kjr = 

j,k j j 

/ \ 1 -a 

E dj+iKj\ 

) 1—ol 

< C"('u;“ 

since O^di < di = w (we cannot estimate this term by W because do = 0) and, because 
of (14. 3D . Y.j>i dj dj+i < ^Y.j>i dj dj < 3W. Finally, we get 

Ej,kMjk < +eT'2||/||2 2 } . (4.9) 

Estimate of the second term in (j4.7p Before following these computations, we need 
to understand what the condition ‘^Bjk H Bji 7 ^ 0” is meaning. Assume that {k,l) G Ij. 
From Qjk Id Qji = 0, we see that either supp pjk or supp is not included in Qjk (resp. 
Qjl). Since, by definition of Ij, we have 9jk > 6ji, this implies that Rjk > dj^\l2. 
Therefore, as we noticed after formula (14.5p . 

Mjk < M'^\\fjk\\\2i^Q.^'^ 

and (14.6p also implies that 

e^,dj+,<e^,dj^,<ce'^-\\fjk\\l.^s^^y 


< 




'j+i 


< C{w + W)^ 
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Now, g3D, the subadditivity of the square root function, the preceding inequality, (031) 
and Cauchy-Schwarz inequality give in turn 


<Ce-^ E \\fit\\l + \\/iM\fi,h 

{k,i)eij 

< I KM\\h(Q) + ,/Eliy‘lli,/Eliy‘lli 

[ k,I ]j k,l 

From Kj < and ||/j'||^ 2 (q) < ll•^llL2(Q)> obtain in the end that 

(4-10) 

j (k,l)&Ij 

End of the proof Let F = e'^^||/|||2- We remind the definition of VF = Wi{f~^, f~). 
One can assume that f~ ^ /+ so that F,W >0. Now, (14.71) . (14.9p . (|4.1Up and the fact 
that J < (7(1 + In re) yield 

M^{u) < <7(1 + I lnu;|) + vj-^F^ , 

where v ■.= 1 — d{l — a) G (0,1) and so a — v = —{d — 1)(1 — a) < 0. Let us fix some 
6 G (0,1) small enough so that 0 < ± J < 1 and — a ± <5 < 0. For some constant c 

depending on (5, one has 1 + | ln?/;| < c(w^ + w~^) and so 

M“(u) < C , 

where the sum is taken over the values of ±1 (+1 or —1) in the right-hand side. Then, 
we make the choice w = W + F^ > 0 for some A = A(a, d) > 0 which will be fixed later. 
Note that all the estimates above are valid only if re < 1. However, if IF -|- F^ > 1 then 
the right-hand side of (|4.1I) is greater than some positive constant and (14.11) easily follows 
from Theorem 13.21 since F[{x) > x. Thus, one can assume that w G (0,1). 

Since 0 < ± <5 < 1, we get and, because —d ± d < 0, 

u — a ± 6 < 0, we have ^ ^ru-a±5 ^-d±5 ^ pX{-d±S) gives 

M^{u) < C |lF''=^^ -k -k IF^'^^^ -(- I . 

We fix A > 0 small enough so that 1 -|- X{—d ± 5) > 0: in this way, all the exponents in 
the preceding formula are positive. Finally, (14.ID follows from the fact that we have IF, 
F < 1 as a consequence of IF, F < -|- F. □ 

, 1 

Remark 4.4. Since minjrc'^ -|- w~^F : tc G (0,1)} = cF <>■+'' and < 1, one cannot 
obtain an estimate of the form M^{u) < C{W -|- F) as expected. However, one could 
improve a bit (14.ip by a better estimate of the number of indices I such that {k,l) G Ij. 


20 













5 A F-convergence result 

Let Q C be a bounded open set and fj, = — fj,~ be a finite measure, where are 

two probability measures compactly supported on fl. We recall the definition of the set 

Aidivi^) = {tt : —>■ : u and V • u are finite measures on tl} 

which is endowed with the topology of weak star convergence on vector measures and their 
divergence. As weak star topology is never metrizable in infinite dimensional Banach 
spaces, the space A4div{^) is not metrizable. Indeed, assume that X is some infinite 
dimensional Banach space such that X' is metrizable. In particular X' admits a countable 
neighborhood basis (Fi)n>i which one can assume to be of the form 

Vn = W '■ l(<y? ;xi)| < en for f = l,...,n} 

for some linearly independent family of vectors (xj)i>i C X and > 0. Then the 
Hahn-Banach Theorem easily provides a sequence (v9n)n>i satisfying (pn{xi) = 0 for all 
z < n G N* and Hv^’nllx' = n. In particular the sequence (<^n)n weakly converges to 0 as 
n —>■ oo which is a contradiction with the fact that {(pn)n is norm unbounded. 

However, every bounded subsets of the dual space of a separable Banach space are 
metrizable for the weak star topology. In particular, for the natural norm = 

|V • + |w|(ri) given by the total variation of u and its divergence, we know that all 

bounded subsets of Aidivi^) nre metrizable: for all M > 0, there exists a metric cIm for 
the weak star convergence of u and V • rt on the set 


AiM{^) = {u £ Aidivi^) : |u|(!d) + |V-u|(!d) <M}. 

In |18| the T-convergence of the functional sequence M“ to M" was proved. Our aim is 
to prove that this property remains true when adding a divergence constraint. Since, for 
u G one has X ■ u £ L?, one cannot prescribe V • tt = /i if /i is not in L^. For this 

reason, we first have to define a regularization of /r. Let {fe)e>0 C be a sequence of 
functions weakly converging to p as measures and satisfying 


[ fe{x)dx = 0 and e^^\\fe\\l 2 —>0 . (5.1) 

JQ 

This choice is going to be useful for the proof of Theorem ll.il For example, we can define 
fe as 

fe ■■= Pe * P, 

where Pe{x) = x) for some compactly supported p £ such that 

f^p = 1 and 7 is still defined as 7 = Now, let us define the functionals 

(resp M°‘) adding a divergence constraint on u G Aidivi^)' 


M^iu) 


M"(n) 

if V • tt = /U, 

+00 

otherwise, 

M-{u) 

ifX ■U = fe, 

+00 

otherwise. 


The main result of this section is Theorem 0 


Theorem. There exists a constant cq such that the functional sequence (Mg)£>o F- 
converges to cqM as e ^ 0. Moreover cq is given by the minimum value for the mini¬ 
mization problem (lOD . 
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We first remind how to build a recovery sequence in the case of a mass 0 flowing on 
a single segment S', i.e. u = To this aim, we need to find a structure close to u 

which is almost optimal for M“. We proceed by a slicing argument: 

Let u be any vector measure in Take some i/ G S^ := {x G : |x| = 1} 

which has to be thought as the tangent vector to S in the case where u = Let 

us consider v = [{u ■ (restriction on of the positive part of u ■ ly) the flux of 

u across the hyperplane = {x G : x • = 0} and assume that f v = 6. Then 

M^{u) can be controlled from below by integrals on subintervals of Mu of the following 
Cahn-Hilliard type energy (see m for physical motivations): 


Ff(u) = e-^i [ u^ + eT'2 f |Vu|2. 

•/ M •/ M 


This kind of models for droplets equilibrium was studied by G. Bouchitte, C. Dubs and 
P. Seppecher in [5] for instance (see also El). Tf(u) can be renormalized through the 
formula: v{x) = ORq'^w^Rq^x), where Re^e = Then, the constraint f v = 6 

turns into f w = 1 and (v) = 9°‘F^{w), where 


F^{w) = I f + f |Vt(;| 

L J M. 


Then, the existence of an optimal profile w is given by 

Lemma 5.1. There exists a profile w G solution of the minimization problem 


mm 


< / u)^ + / Irc'p : w & and I re = 1 1 . (5.2) 

VM M M J 


Moreover, w is compactly supported, Lipschitz continuous on M and C°° inside its support, 
i.e. on the open set {w > 0}. Last of all, it is possible to choose w such that it is even 
and non-increasing on M'*'. 

Remark 5.2. Although we restrict to the two dimensional case, some works by G. Bou¬ 
chitte, C. Dubs and P. Seppecher (see |13| 1 suggest that Lemma l5.ll and Theorem 11.11 
could be generalized in every dimension as well. However, the aim of this paper is to use 
the tools of section U] so as to establish the P — lim sup property for functionals (with 
divergence constraint) and, from the point of view of the complexity of the proof, this 
is independent of the dimension. Since the P — lim inf property was only established in 
the 2D case in [18| . we prefer to stay in this framework. Actually, the difficulty to prove 
a P-convergence result of (resp. M“) to M (resp. M“) in every dimension would 
concern the P — lim inf part and this is not the purpose of this paper. 

Remark 5.3. Note that the minimum value in (15.21) is related to the best constant in the 

1 —/3 3 

one-dimensional Gagliardo-Nirenberg inequality /jg |w| < C 


C 


= inf 



2+/3 


u 




loc\ 


and 


/ 

Jr 


u = 1 


Proof. First notice that there exists a finite energy configuration, i.e. w G 1/^(M,M) such 
that Fl^{w) < -|-oo. Indeed, every compactly supported and nonnegative function has 
finite energy. Let take a minimizing sequence, i.e. {wn)n C H^{M) C C^(M) such that 
F^{wn) cg. One can assume that Wn is even and non-increasing on M"*". Indeed if w* 
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stands for the spherical rearrangement of some w G then one has \w*\^ = 

and the classical Polya-Szego Theorem states that the spherical rearrangement 
reduces the Dirichlet energy of w: 


IK)?< [ 

Jr 


..'1 2 


W 


In particular F^{w*) < F^{w) as announced. Since {w!^)n is bounded in one can 

assume that it weakly converges in L^(M). Moreover, as Wn > 0 a.e. and = 1, 

{wn)n is bounded in Thanks to the Poincare-Wirtinger inequality, one deduces 

that {wn)n is bounded in Up to extraction, one can assume that {wn)n weakly 

converges in Let call w G the limit. In particular, {wn)n strongly 


converges to w in and so w is even, nonnegative and non-increasing on M"*". 

Moreover, the Fatou lemma and the weak convergence of w'^ yields 

F^{w) < liminf F^{wn) = cr . 

In order to prove that re is a global minimizer it remains to prove that w satisfies the 
constraint f w = 1. Indeed, from the Fatou Lemma we can only deduce that f w < 1. 
One has to prove the strong convergence of Wn in L^(M). Since Wn converges in it 
is enough to prove that the sequence (wn) is tight. Let R > 0. For all n > 1, since Wn 
is non increasing on [0, i?], one has ||rrJn||L°°({a;: | 3 ;|>J?) = WniR) and Markov’s inequality 
yields Wn{R) < ^ = ^- Hence 


f Wnix) dx < Wn{R)^ ^ f W^{x) dx < {2R)^ ^F^{Wr, 
J\x\>R JM. 


< 


c 




for some constant C > 0 non depending on n which implies that {wn)n is tight since 
1 — P > 0. Now, let check the regularity of w: Lipschitz continuous and smooth inside its 
support. Note that we already know that w G thanks to the Sobolev embedding 

F[i^^(R) C Let check the regularity of w inside its support. Since w is even and 

non-increasing on M"*", the set {x : w{x) > 0} is an interval {—R,R) for some R > 0. 
Since tc is a minimizer of the minimizing problem (15.2p . w also satisfies the following 
Euler-Lagrange equation 


Vx G {—R,R), —w''{x) + j3'w{xY ^ = A 


(5.3) 


where A G M is the Lagrange multiplier associated to the volume constraint f w = 1. 
Note that A = —w"(0) -|- > 0. Indeed, since x = 0 is a global maximum of w, 

w'(0) = 0 and w"(0) < 0. From (|5.3p . one deduces that w is smooth on (—R,R). Now, 
multiplying (|5.3p by w' and integrating it on [0, x] or [x, 0] yields 

/ / A 2 

Vx G (—R, R), —— = w{x)^ — \w{x) + Ar(;(0) — . 

Since w is bounded on M, we deduce that w is Lipschitz continuous on M. Last of all, 
we prove that w is compactly supported, i.e. R < oo. Assume by contradiction that 
R = oo. Then / G / > 0 on M and integrating (15.31) yields 


Vx G M, w'{x) = f {/3w{y)^ ^ — X)dy 

Jo 


10 


Since w{x) —> 0 and /3 — 1 <0, w{x)^ ^ —> -|-oo and the right hand side of the 

X^OO X^OO 

preceding equation goes to -|-oo as well. Thus w'{x) —> -|-oo which is a contradiction. 

X—)-CO 

□ 
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Now let wo^e be defined by W 0 ^e{x) = 9Rq'^w{RqIx) where w is the optimal profile 
of Lemma 15.11 (satisfying all claimed regularity and symmetry properties) and let us 
introduce the kernel associated to we^e, given by the following lemma 

Lemma 5.4. There exists a bounded and compactly supported radial kernel M"*") 

such that is the projection of pg^^ on the axis (xi = O).' 


Ujpg^^{x) dx = Wg^e{x 2 ) dX 2 , 


where 11^ stands for the projection on the second variable, dx (resp. dx 2 ) is the Lebesgue 
measure on (resp. and nj/i stands for the pushforward of some measure p by 
n : —>■ M. Moreover, one can choose pg^^ of the form pg^e{x) = Rf'^p{RQ^x) for some 
pGL“(M,M+). 


Proof. After the renormalization pg^^(x) = ^x), it remains to hnd p satisfying 

njp(x)dx = w{x 2 ) dx 2 . It is not very difficult to see that a radial solution is given by 
the formula 



(5.4) 


Details are left to the reader. Let justify how (15.4p implies that p is bounded. Since w is 
compactly supported, there exists i? > 0 such that w{x) = 0 for |x| > R. In particular, 
p is compactly supported on B{d, R). Then, since t(;'(0) = 0 and since w' is bounded on 
M and smooth around 0, one has |'u;'(x)| < C\x\ for all x G M and some C > 0. Hence, 
there exists a constant C > 0 such that for all x G such that |x| =: r G [0, R), 


p{x) 



—w'{r.s) ds 

— 1 



rs ds 

7P^ 


< c 




s ds I 


which is bounded since s —?■ ^ is integrable on [1, 2] and bounded on [2, +oo). □ 

As a consequence, in the case where u = OTL^^, a recovery sequence, i.e. a sequence 
(ug) such that Ug ^ u in Aidivi^) and M°‘{u) as e —>■ 0, is obtained as 


Ue = P9,e * U. 

In the case of a hnite energy configuration, i.e. u G Mdivi^) such that M°‘{u) < oo, 
thanks to classical properties in the theory of L-convergence, it is enough to find a 
recovery sequence for u belonging to a class of measures which are dense in energy. 
Thanks to the work of Q. Xia in |23) (see also |18|L we know that the class of vector 
measures concentrated on hnite graphs is dense in energy so that one can restrict to this 
case. In |23| . the branched transportation energy was in fact dehned by relaxation of its 
restriction to the set of vector measures concentrated on a graph. The rectihability of 
hnite energy conhgurations and the Eulerian representation (|2.1I) was discovered later in 
[24] (see also El). This density property was used in [18] to prove the T-convergence of 
M“ toward M". In the setting of functionals with divergence constraint, we need the 
following lemma: 

Lemma 5.5. Let u G Aidivi^) be such that M°^{u) < oo. For all A > 7 , there exists a 
sequence (u^) C Hq{Q) converging to u in Aidiv(^} such that 

—> coM°‘{u) and e'^||V • tt£|| 2,2 is bounded. 

£^0 
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Before proving this statement, we are going to investigate the case where u is concen¬ 
trated on a finite graph. First of all, let us give some details on what “a vector measure 
concentrated on a finite graph G” is. Let G = {V{G),E[G),9) be a weighted directed 
graph: V{G) C O is a finite set of vertices, E{G) is the finite set of oriented edges 
e = (e, Te), where e = [oe, 6 e] C and Tg is a unit vector representing the direction of e, 
and 6 : E{G) —(0, -|-oo) is the weight function. Then the “vector measure associated to 
G” is given by 

UG= ^ 0{e)Te d^|g . 
e={e,Te)&E(g) 

These measures uq belong to Aidivi^), he. V • uq is a measure, and they are called 
“transport paths” (see Definition 2.1 in |23|i. When n is a transport path, we have the 
following lemma: 

Lemma 5.6. Let u = uq £ Aidivi^) for some weighted directed graph G. Then, there 
exists a sequence {us)e>o converging to u in Aidiv{^) ond a constant G depending on u 
such that, for e small enough, G and 

L Jq \ue\ < + Ge"^, 

In 

3. ||V • rt£||x,2 < G, 

I \M^iue)-coM^{u)\<Ge'r. 

Proof. By definition, such a vector measure u can be written as a finite sum of measures 
Ui = 9i Ti concentrated on a segment Si C 0, directed by r, with multiplicity 6 i 
for i = 1,...,/. We first define a regularized vector fied Ve by := where 

Vi = P 9 i,£ * Ui- Then, for e small enough, is compactly supported on D and satisfies 

f /n 

\\M^{ve)-coM^{u)\<GeA 

The first statement is a consequence of the fact that j = 1 and the inequality 
11/ * pWl^ < II/IIli Ia*I(^) foi' / G Cc{Ll) and for a finite measure p on D. For the 
second statement, by definition of the kernel pg^^ we know that, out of the nodes set 
N = UiSupp(V • Vi), 

M“(u„iV'=) = coM“(u,A^‘=). 

As a result, we just have to estimate these energies on N which is a finite union of balls: 
the supports of pg^^^ recentered at each end-point of the segment Si. Since the radius of 
these balls is of the order of e'^, this immediately gives the fact that M^[u,N) < Ce'^ 
for some constant G > 0 depending on u. For the sake of simplicity, in the rest of this 
proof, G > 0 will denote some constant depending on u which is large enough so that all 
the inequalities below are true. We are going to prove that 

M“(m, N) + M^{ve, N) < Ge^. 

It remains to estimate M^{vs, N). Since M^{ve,N) < I Yli ^ei^i, N), it is enough to 
estimate M^{vi,N). But ||'yj||Loo(^) = ||/06»i,£ *'Wj|| 2 ,oo(jv) < C£~‘^"i\\p\\L^\ui\{Ni), where 
Ni\= N + suppp 0 .,£ ■.= {x + y : x e N, y e suppp^i.e}- Note that suppp 6 )i,e is a ball 
centered at the origin with radius smaller than Ge'^ so that W is a finite union of balls 
with radii smaller than Ge'^ as well. In particular, using the fact that Ui = OiTiTL^^,, 
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we get \ui\{Ni) < Ce^ and so ||nj||x,oo( 7 v) < Ce . Similarly, one has ||Vni||ioo(^) = 
II< Ce~'^'^\\S/p\\L'=°\ui\{Ni) < Now, the definition (|2.5I) gives 

M“(n„lV) = e“+^ [ |Vn,|2 + e“-i [ \vef < \N\{e^+^\\Vv,\\lo. + 

J N J N 

From the inequality |A^| < Ce^'^ and the equalities a = 87 — 1, /?7 = =47 — 2 , 

we deduce 

M^{Ve,N) < Ce2T'{e37-47 +^37-2-/37} < 2Ce^ 

as required. In order to construct an approximating vector field with controlled diver¬ 
gence, we need to consider := where G Hq{N) is constructed as follows: 

The node set N, defined above, is a finite union N = UJ^=i where each node Bj 
is a ball centered at the end-point ai of some segment Si = [ai,bi]. Let assume that £ is 
small enough so that these balls are non-overlapping. Then, on each node Bj, gj := 
is a finite superposition of kernels like recentered at Cj, the center of Bj. In particular 
\\9j\\L^{Bj) ^ C£~'^ and pj = V • Ug = (V • u){Bj) =: 9j. If Oj = 0, then Theorem 
allows to find Wj G H^^Bj) satisfying ^ -Wj = pj and \\wj\\fji(^^.-^ < Ce ^. If 9j 7 ^ 0, 
let say 9j > 0, we rewrite pj as pj = p^ — g~ = Xp^ + (1 ~ X)g~^ — g~ where g^ (resp. g~) 
stands for the positive part (resp. negative part) of g and A G (0,1] is chosen such that 
(1 “ = !b9~^ i-®- ^ Ibj 9+- Applying Theorem [331 we get Wj G H^iBj) 

satisfying V ■ Wj = {1 — \)g'^ — g~ and || < C £~'^. Let us define We = J2j 

and Ue := — w^. Since |V • Ue\ = /g. \gj — V • rcjl = A f^. g^ = 9j for all j, we have 

f \V-Ue\= [ |V-Ue| < = |V-u|(fI). 

Jn J N , 


Moreover, to estimate ||V • 11^,2, note that ||V • rt;£||g2 < ||t(;e||gi < C£~'^ and, because 

V • Ue is only composed of a finite sum of translated kernels of the form pe^^£, ||V • r;e||g2 < 
C£~'^ as well. In particular e"*'||V -^£11^2 is bounded. Then, the Poincare inequality yields 


W£ 


Il2 — '^\\wj\\L2{Bj) < C'X^e’^l|V^illi,2(B^.) < C 


since Bj is a ball of radius C'e'^. Consequently, by the Cauchy-Schwarz inequality, we get 

[ {Uel < [ l^^el + [ Iwel < |n|(n)|iV|^/^||r(;£||g2 < |n|(n)C'e'’'. 

Jn Jn Jn 

Similarly, by a Holder inequality, we have 

\we\^ < liVl^lIrn^llJ^ < in£^'^)\^\\w£\\l, < 


J N 

Once again, since a = 87 — 1 and jJ'y = 47 — 2, we deduce 

M“(rc£) = e“+i [ iVtCsp + e“-i f + e«-i+7(2-/3)} ^ 2C£^. 

Jn Jn 

Since M°‘{u,N) < Ce'^, we get M^{us,N) < 2 [M^{v£, N) -|- M^{we, N)\ < Ce^ which 
finally gives 


\M^{ue) - coM“(n)| = \M^{u£,N) - coM^{u,N)\ < Ce\ □ 
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Proof of Lemma \5. ,5i First fix u E Aidiv{^) and construct a sequence {un)n>i converging 
to u such that Un = is a vector measure associated to some weighted directed graph 
Gn C Ll and M°‘{un) converges to M°‘{u). Since (u„) weakly converges in A4div(^}, the 
total variations of both measures Un and V • Un are bounded by some constant M > 0. 
In the following, we use a metric d on the space Extracting a subsequence if 

necessary, one can suppose that the two following estimates hold 

d{un,u) < 2-^-1 and |M"(un) - M“(n)| < 2 -"-^. 

For each n > 1, let Ue^n be a sequence converging to as e —)• 0 and satisfying all 
properties in Lemma [5. 6 1 for some constant C = Cn- Then, one can construct by induction 
a decreasing sequence {en)n>i 0 such that for all n > 1 and s < Sn, and 

1- ^e,n € 

2 . d{Ue,n,Un) < 2 “"“^ 

3. \M»{ue,n) - coM“K)| < 2—\ 

4. s^~'^Cn < 1 SO that e'*' ||V • u^^nWi^ < 1- 

Indeed, assume that > 0 satisfies all the asked properties. Then, one can find £n+i G 
( 0 ,e,i) small enough so that 

* < 2 “”“^ thus implying the first and third properties (see properties 1 ., 
2. and 4. in Lemma [5.61) . 

* Cn+i 1 which is possible since A > 7 

* and d{ue^n+i,Un+i) < 2“"'“^ for all e E (0,en+i) which is possible since Ue^n+i 
converges to Un+i in {M.M+i{^)id) as e ^ 0. 


Now it is quite straightforward that the sequence {us)e>o defined by 


Ue 


Ue,l if £ > £ 1 , 

Ue,n if Gn+i < £ < Gn for some n > 1 , 


satisfies all the properties of Lemma 15.51 


□ 


_ p 

Proof of Theorem M . 1[ It is already shown in [18| that M“ —> coM°‘. We just have to 

prove that the F — limsup property still holds when we add the divergence constraint. In 

other words, it remains to prove that for all u E Mdivi^) such that V-rt = /a, there exists 

a sequence (u£)£>o C M.div{fd) weakly converging to u as measures, satisfying V -Ug = /^ 

(so that (ug) also converges in Aidivi^)) and M“(ue) —> cqM^{u). 

e—>-0 


First of all, take a sequence (u£)e>o C P[q{LI) converging to u given by Lemma 15.51 
for A = ^: M^{us) —>■ M^{u) with e'’'||V • Us\\l 2 bounded. Then define := /e — V • u,. 
the residual divergence. Note that = 0. Indeed /e = 0 by assumption (see (15.IM 
and f^V -Us = 0 since E Lfg(n). Moreover, our assumptions on the sequences (u£)£>o 
and (fe)£>o yield 


S —^0 


Indeed, we know that the same estimate is satisfied by (/£)£>o thanks to (15.11) . Moreover, 

K "y 5'7 "7 

since 87 = 72 and A = one has ||V-M£||i 2 = £ 2 e~^ ||V-M£||i 2 < Cez — 5 - 0. Moreover, 
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since and V • weakly converge to /r as e goes to 0, we know that weakly converges 
to 0. Let g'^ (resp. gj) denote the positive (resp. negative) part of g^. In order to satisfy 
the divergence constraint, we may correct with a vector field w^, given by Theorem 
14.11 (together with Remark I4.2p , such that V ■ Ws = ge, 

M^iwe) < H ^ 0 (5.5) 

and \\we\\Mdiv(^) bounded, where H{x) = C{x + x^) for some C > 0 and 6 G (0,1). 
We deduce that (ws) is relatively compact in Aiciivi^)- From (15.5p and the F — liminf 
property, this implies that converges to 0 in Aidivi^)- Now, by construction, = 
Ue + We satisfies V • = /g, —)• tt in A4div(^) and M^(ve) —coM"(u). Indeed, this 
last limit is a consequence of □ 

Lemma 5.7. Let Q be some bounded open set in d > 1. Let (ue), (ug) C be 

two sequences such that M^{ue — Ve) —)• 0 and assume that M“(ug) is bounded. Then, 

|Mf(ug)-M“(ug)|^ 0 . 

Proof. Let > 0 be some constant. For all real matrices A and B of size d x d, by the 
Young inequality, we have 

|y4 + R|2 = |A|2 + |B|2 + 2A: B <{1 + ufAl"^ + (1 + l/i^)|^P- 

Writing Ue = Ve + Ue — Ve, we use the preceding inequality for A = Vve, B = V(ug — Ug) 
and the subadditivity of x ^ |x|^ to get 

Mf(ug) = f \uef + [ \Vue\^ < (1 + i^)M“(ug) + (1 + l/i^)M“(ug - Ug). 

Ju Jo, 

Since M“(ug) < C for some constant C < +oo, we deduce that 

M“(ug) - M“(ug) < Cz. + (1 + l/z/)M“(ng - Ug). 

For any value of e such that Ug 7 ^ Ug, let take v = \/M^{ue — Ug) > 0. Hence, taking the 
lim sup when e —)• 0 , one gets 

hmsup{M"(ttg) — M^{ve)} < C'lim sup Mf‘{ue — Ug) = 0. 

£—)‘0 £—^0 

Since M"(ug) < 2[M“(ug)+M“(ug—Ug)] and M“(ug) is bounded, we deduce that M°{ue) 
is bounded as well. Then we can apply all the preceding computations exchanging ttg 

and Ve to get hmsup{M“(ug) — M"(ug)} < 0 which concludes the proof. □ 

£— ^0 
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